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abstract 

Density fluctuations of fluids with negative pressure exhibit decreasing time behaviour 
in the long wavelength limit, but are strongly unstable in the small wavelength limit when 
a hydrodynamical approach is used. On the other hand, the corresponding gravitational 
waves are well behaved. We verify that the instabilities present in density fluctuations are 
due essentially to the hydrodynamical representation; if we turn to a field representation 
that lead to the same background behaviour, the instabilities are no more present. In 
the long wavelength limit, both approachs give the same results. We show also that this 
inequivalence between background and perturbative level is a feature of negative pressure 
fluid. When the fluid has positive pressure, the hydrodynamical representation leads to 
the same behaviour as the field representation both at the background and perturbative 
levels. 

PACS number: 98.80.Bp, 95.35. +d 

keywords: cosmology, perturbations theory, large scale structure. 



1 Introduction 

Fluids with negative pressures have become of utmost importance in modern cosmology. 
They were first considered in the context of the inflationary scenario in which the early 
universe has a very short period of accelerating expansion |JJ. This inflationary phase 
solves many problems of the standard model which are connected to specific choices of 
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initial conditions, like the flatness and horizon problem. At the same time, the origin 
of the seeds of the large scale structures observed today in the universe has a natural 
explanation in the inflationary scenario, that considers them as quantum fluctuations of 
a scalar field in a de Sitter background. The inflationary phase must end up with a 
transition to the radiative phase of the standard model. In many aspects, the inflationary 
scenario is plagued with problems of arbitrary choice of parameters of the microphysics, 
but its great success with the above mentioned problems led the community to consider 
it as part of the standard scenario. 

More recently, the observation of the supernova of high redshift led to the preliminary 
conclusion that the Universe today is in an accelerating expansion [0, [3] . This is a surpris- 
ing result since there was little doubts that the Universe was desacelerating and that the 
desaceleration parameter q was positive. If the results of the supernova observations are 
confirmed, the energy of the Universe today must be dominated by a fluid with negative 
pressure. It could be a cosmological constant, but other cases are not excluded, like a 
fluid of cosmic string or domain walls 0, |5|, [| or some scalar field, with a peculiar kind 
of potential, called quintessence [171 Bl [| [II], [TT[] . 

In many pratical cases, these fluids with negative pressure may be represented by a 
perfect fluid with a barotropic equation of state p = ap, with a < 0. To be more precise, 
in order to have an accelerating universe the strong energy condition must be violated 
and a < — |. A fluid of cosmic string leads to a = — |, representing the limiting case 
between an accelerating and desaccelerating universe; a — — | and a = — 1 represent 
respectively a fluid of domain wall and a cosmological constant. The quintessence fluid 
only approximativelly can be expressed by a barotropic fluid with a < 0. 

In [J2j perturbations in fluids with negative pressure were studied in the hydrody- 
namical representation. It was found that when the strong energy condition is violated, 
there are instabilities in the small wavelength limit. This result can be understood by 
remembering that in this limit the expansion of the universe plays no important role: the 
negative pressure at this limit acts in the same way as gravity and nothing can prevent 
the collapse. This can be easily seen in the case of a cosmic string fluid for which the 

1±a/ i+— 

density contrast behaves as A oc t V 3 . When n — > oo divergences appear, leading to 
the instability of the background model. In the cosmological perturbation theory, we are 
generally interested on the unstable modes, but that are not divergent. From here on, 
instabilities will refer to these undesirable divergent modes. 

This study of density perturbations for a fluid of cosmic string has been extended 
in [IB], where a two-fluid model was considered, one of them being the string fluid and 
the other ordinary matter with positive pressure. Special attention was payed to the 
case of a closed spatial section. From the point of view of the background, such models 
have many interesting features, mainly conected to the horizon problem and to the age 
of the universe. However, the study of scalar fluctuations around this background shows 
the presence of instability in the small scale limit, as in the case of the one fluid model. 
On the other hand, the study of gravitational waves for these models |14| reveals a very 
regular behaviour because gravitational waves are mainly connected with the scale factor 
behaviour, being quite insensitive to the matter representation. This indicates that the 
instabilities detected in |l2|, |l3| could be due to the hydrodynamical representation, which 
could disappear in a more fundamental approach. 

In fact, the hydrodynamic representation for fluids with negative pressure is frequently 
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a very poor approximation. Negative pressures appear in situations involving phase tran- 
sitions in a primordial universe (topological defects) or a fundamental self-interacting 
scalar field. The exact formulation involves consequently fields, and a representation us- 
ing a fluid with a barotropic equation of state only in very special situations may be 
employed. The employment of a perfect fluid representation, mainly when fluids of neg- 
ative pressure are involved, can be viewed as a practical simplification which in many 
situations gives the same results as those that could be obtained by employing a more 
fundamental field representation. 

The main point is that the equivalence of a hydrodynamical representation with a 
more fundamental one is very restricted and it can lead to complete misleading results 
depending on the problem treated. The instability in the small wavelength limit quoted 
above is an example. Fluids with negative pressure should have a field representation 
where their main features could be retained. The representation we will investigate here 
involves the more reasonable coupling of gravity with a scalar field with a self interacting 
term. 

The aim of this paper is to show that the scalar field representation can more conve- 
niently retain the features we could expect from fluids with negative pressure, mainly for 
those that are interesting for cosmology, as the objects resulting from phase transitions, 
like cosmic string. We will concern mainly with the consequences of these different repre- 
sentations for an analysis of perturbations around a homogenous and isotropic expanding 
universe. We verify that when we use a field representation, with the same behaviour 
for the scale factor as the corresponding hydrodynamical representation, the instabilities 
present in the latter case are absent in the former one. This is due essentialy to the fact 
that when we pass from a hydrodynamical representation to a field representation of the 
fluid we also pass from an Euler's type equation to a Klein-Gordon's type equation, and 
there is no correspondence when the pressure is negative. On the other hand, for large 
scale perturbations, both approachs lead to the same results. 

We begin by analysing a stiff matter fluid which can mimic a scalar field in a very 
simple way. As it is well known, a free scalar field minimally coupled to gravity reproduces 
the stiff matter equation of state. We will show in section 2 that for the stiff matter 
and free scalar field models, the agreement between the results occurs not only at the 
background level but also at the perturbative level. In section 3 we extend this analysis 
to a perfect fluid with arbitrary equation of state p = ap: we determine the corresponding 
field representation, showing that, at perturbative level, the equivelence between these two 
approachs remains only when a > 0. In section 4 we review briefly the results obtained in 
fT3] , |i~4f and we discuss the possibility of a field representation for this two fluid model. In 
section 5, we review a model of variable cosmo logical constant that leads, from the point 
of view of the evolution of the scalar factor, to the same behaviour as a cosmic string 
fluid. In section 6, we perform a perturbative analysis of the variable constant model, 
and show explicitly that they are free of instabilities, both for scalar perturbations and 
tensor perturbations. In section 7, we discuss the results. 
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2 Free scalar field model 



The most simple gravity model with a scalar field is described by the action 

s= [ V^g[R~^ a( p, a }d 4 x . (i; 



It represents a free scalar field minimally coupled to gravity. The field equations obtained 
in accordance with the principle of least action by varying S with respect to the dynamical 
variables are: ^ ^ 

Rnv - -g^R = 4>-,A\v ~ ^9^-A P , ( 2 ) 
□0 = . (3) 

In a Friedmann- Robert son- Walker flat space-time (FRW), the metric describing the four 
dimensional geometry is 

ds 2 = dt 2 - a 2 (t) f^dx 2 + dy 2 + dz 2 ^j , (4) 
and the field equations take the form 

'(;)" - • < 5 » 

a f a\ 2 d) 2 
2-+ - = -\ , 6 

+ 3-0 = , (7) 
a 

where the overdot denotes the derivative with respect to the time coordinate t. The 
equations (HJ^J^) are not independent due to the Bianchi identities. 

On the other hand, in an universe filled with a perfect fluid we have the field equations 



Rnu — t;9livR — SttGT^ 



1 

T 

T£ = , (9) 

where 

Tp„ = (p + p)u^u v - pg^ , (10) 

with p = ap, a being a constant. The most common values of a of cosmological interest 
are (pressurelless matter), | (radiation) and 1 (stiff matter). Negative values of a 
has acquired increasing importance due to the inflationary paradigm and the new results 
coming from the supernova type la observations, as it was discussed before. Topological 
defects also require a negative equation of state. Again, equations (§|§) are connected by 
the Bianchi identities. 

In the same FRW flat background, (||,[|) take the form 

= 87rGp , (11) 

a ( d\ 2 , , 

2-+ - = -8nGp , 12 
a \aj 

p + 3(1 + a)-p = . (13) 
a 
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It is a straightforward to check that the equations (|||6|fl) an d equations (|IT|, |i2| , |i"3| ) 
permit the identification 

i.e., the "scalar field fluid" behaves like a stiff matter in the hydrodynamical approach. 
The scalar field sound velocity, in this case, is A = p^/p^ = 1. The scale factor both in 
the free scalar field and hydrodynamical stiff matter cases behaves as a oc t 1 / 3 . 

The evaluation of the perturbed quantities follows the well known approach of Lifshitz 
and Khalatnikov P~5l [TBI. It can be treated either with the synchronous gauge or the 



gauge-invariant formalism. Here we choose to work in the synchronous gauge formalism, 
where h^ = 0. 

We study first density perturbations and then gravitational waves for this free field 
model. 



2.1 Density perturbations 

Introducing small perturbations around the background solutions, the perturbed equa- 
tions for the scalar-tensor model read 

h + 2-h = 4<j)5<j) , (15) 
a 

(I • Tl 2 1 • • 

5<j> + 3-50 + — 5<j) - -h(j> = , (16) 
a a 1 I 

where h = hkk/a 2 and n 2 comes from the Helmholtz equation V 2 Q + n 2 Q = 0: the scalar 
functions Q(x k ) are the eingef unctions of the three-dimensional Laplacian operator. 

In order, to solve equations (|15| , |T6|) it is more convenient to work in the conformal time, 
dt = adrj. The scale factor behaves as a oc r/ 1 ^ 2 . In terms of this new time parameter the 
solution of the perturbed equations is given by |IT 

S(p = V~ 3/2 J 77 3/2 (ci(n)Ji(n??) + c 2 {n)N 1 {nr])\dr] , (17) 

where J\ and N\ are, respectively, the Bessel and Neumann functions of the first order, 
and Ci(n), c 2 (n) are two arbitrary constants. 

We need to verify if the evolution of the perturbations are specified by the equation 
dpc/, = a5pj> as it happens with the background evolution showed previously. This should 
imply that the adiabatic approximation is verified here. To do this, we consider the 
perturbation of the equation (0) and fllUD 

+ ^h^cf>, p <f>, a —g^g^iS^p <f>, a +0 )/3 6<b, a ) , (18) 

5T^ = (8 P4> + Sp^WU" + (p + p^Wir + U»5U V ) - Sp+gT + P^ v . (19) 
Using the synchronous gauge condition, we have 

5G 00 = <p5<P , (20) 
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SGij = --h i:i ft + a 2 6 i:j <p6<p , (21) 

5T 00 = 5p^ , (22) 

STij = -p<phij + cfdijdptp . (23) 

If we consider the perturbed Einstein equations 5G^ V = 8tt GST^ and the values of p^ 
and obtained from equations ( TTTp and QT^), we have that 

= . (24) 

In the hydrodynamical approach, the solution for the Einstein perturbed equations 



with p = p and 5p = Sp, leads to the expression [fL2]| 

A = — = rr 3/2 | ^ B/2 (di(n)J (rM7) + d 2 {n)N {nr 1 )}dr 1 , (25) 



where rfi(n) and 0^2(71) are the integration constants. Remembering that = f-y and 
bp^ = 2^r~, the quantity A<^ = computed from ( [TTD reduces to (p5|), using simple 
recurrence relations for Bessel's functions. 

Hence, in this simple model where the matter is described by a scalar field, the back- 
ground and perturbed equations can be put in a barotropic form. The "velocity of sound" 
is the same as the one defined by A —p^jp^ in agreement with fi~8f , where it is 



shown that this result corresponds to the low-frequency regime of the scalar field oscilla- 
tions. 

2.2 Gravitational waves 

Here, the perturbed equation of gravitational waves is: 



h " _ 2-ti + 



a" a' 2 
n 2 -2— + 2— 

a or 



h = , (26) 



where hij = h(i])Qij, being a tracelless transverse eigenfunction on the spatial section, 
and primes denote derivative with respect to conformal time drj = adt. 

After inserting the backgroung solutions, we obtain the solution of the equation (|2"6"|). 
as follow: 

h = r}\ei{n) J\{nrj) + e 2 Ni(riT])j , (27) 

where e\(n) and e 2 (n) are the integration constants. This solution is valid for both 
represent at ions . 

It is easy to verify that the above solution is well-behaved and stable. 

3 Field representation for an arbitrary barotropic equa- 
tion of state 

Let us return to the Einstein's equations coupled to a perfect fluid, with a barotropic 
equation of state p = a p. Solving the Einstein's equation for a flat spatial section, we 
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obtain for the scale factor a = aot 3(1+a) . Let us now consider a self interacting scalar field 
coupled to gravity. The Lagrangian has the form, 



L 



-9 



R - <j>, p <jyP + 2V(<i>) 



(2* 



where the potential V(<f>) represents the self interaction term. The field equations are 



R 



1 



g^R = <f>.p<f>. v - -g^4> ]P (p' ,p + 9nvV{4>) 



□ 



= - W) 



(29) 



where means partial derivative of the potential with respect to 0. 

This scalar-tensor model may lead to the same behaviour of the scale factor as in the 
perfect fluid model provided that the potential takes the form, 



v{<(>) 



2 (1 -a) 



exp(W3(l + a)4>) . 



3(1 + a) 2 

Indeed, for a FRW background, this potential leads to the solutions 

2 



a = a t 3(1+a) 



3(1 + a) 



hit 



(30) 



(31) 



For the case of a cosmic string fluid, a oc t, the potential is V((p) = —2e~^^. 

We now turn to the perturbative level. The perturbation of the scalar-tensor model 
leads to 



h" + -ti 
a 



+ 2 



2-6' 



b" + 2-6'+[n 2 + V H a< 



-S'ti 
2 y 



(32) 
(33) 



where we follow the same definitions as before and we have employed the conformal time. 



Combining equations (|32| , |33|) , and inserting the background solutions, we obtain the 
third order differential equation 



+ 



7 + 3a i 



+ 3 



1 + 3a 7] 
1 + a n 2 



2 - 24a - 18a 2 



1 + 3a) V 
2 1 



a 



1 + 3a rj 



1 + 3a) 2 rf 



(34) 



One solution for this equation is 5(p = r/~ n- 3 «, which is given by the residual coordinate 
transformation freedom characteristic of the synchronous coordinate condition. Using this 
known solution and some suitable transformation, we can reduce the third order equation 
to a second order Bessel equation. The final solution for 5<j) is 



n 



c±i]2 J ±v {nrj)drj 



(35) 



where c± are integration constants, that in general depend on n, and J± u is a Bessel 
function of order v = f + , 3 „ a , . 

2(l+3«) 



The solution for the density contrast in the case of the hydrodynamical representation, 
with p = ap is well known [p^j . It can be written as 



A = i] 3 i+3a J rjUc+J^y/anr]) + c-J-^y/anrf^di] ,a>0 , (36) 
A = V -^J V i(c + I^n,) + c_K^n,))d V , „<0 , (37) 

where now 1^ and are the modified Bessel functions, and p — | fiiff^ • I n or der to 

connect both results we must remember that = where p<j> = \ -\- V((f>). Using 
some Bessel's recurrence relations, as in the previous section we find 

A^ = C±?T 3 ^ j(l - -)v^J±n(nv) + J V*J±i*{™j)dv\ ■ ( 38 ) 



This expression differs from (p6|) by the factor a in the argument of the Bessel function, 
and by the first term in fl38|). However, when a > solutions (|36| , |38| ) have the same 
behaviour in the large and small wavelength limit. For a — 1, the equivalence between the 
two approaches is complete, as can be easily checked by comparing (|38| ) with (^) for this 
special case. For a < 0, the correspondence between solutions ([37],[3|) occurs only for n — > 
(large scale perturbations); for n —>■ oo (small scale perturbations) the hydrodynamical 
representation exhibits strong instabilities while the scalar field representation displays 
accoustic oscillation. 

Finally we remark that both representations give the same behaviour for gravitational 
waves. 



4 Perturbations in cosmic string fluid 

An important particular case of the scalar-tensor model develloped in the previous section 
concerns the case of cosmic string. A cosmic string fluid may mimic a perfect fluid with 
an equation of state p — — |. In [[HJ the fate of density perturbation in fluids with 
negative pressure has been studied. The main conclusion was that, in the long wavelength 
approximation, there are only decreasing modes when the equation of state is such that 
the strong energy condition is violated; for small wavelengths, instabilities can arise due 
to the fact that the pressure contributes in the same sense as gravity. For the equation 
of state of a cosmic string fluid displayed above, density perturbations behave as 

A = r !±\A+? . (39) 

Hence, in the small wavelength limit, n — > oo, instabilities come to scene. 

A scalar-tensor representation of the cosmic string fluid is given by a scalar tensor 
model with V((f)) = 2 exp(±v / 20). A perturbation analysis as it was performed in section 
3 leads to the following expression for the perturbation in the scalar field (the integration 
procedure follows very closely to that sketched in the previous section): 

A^ oc r 1±VT ^ . (40) 
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In the large wavelength limit both representations give the same results; in the small 
wavelength limit the hydrodynamical representation display instabilities while the scalar- 
tensor model exhibits accoustic oscillations. 



In |13| it was considered the case of two non interacting fluid, one of them represented 
by an equation of state like p = — (1/3) p while the second fluid is the ordinary matter 
with a barotropic equation of state p = a p. This may represent in more realistic way the 
presence of fluids with negative pressure in the universe. The equations of motion are 



3(-) 2 + 3^ = 8nG( Pm + p s ) , (41) 
a a z 



a ,a, 9 k 8irG , „ s , in > 

2- + {-) 2 + - = —{ Ps -3ap m ) , (42) 
a a or 6 

p m + 3-(l + a) Pm = , (43) 
a 

p s + 2-p s = . (44) 
a 

In this equations p m and p s denote the ordinary fluid and stringlike fluid densities 
respectivelly. 

The background solutions are: 

a = a sinh 2 (^-?7) , when a = ; (45) 

a = a sinh( v /7?7) ,when a = 1/3 ; (46) 



a = aoysinh( 2-^/777) ,when a = 1 , (47) 
where 7 = \^nGp 0s — k\, p 0s is defined as the ratio p s = rj being the conformal time. 

4.1 Density perturbations 

Perturbing the Einstein's equations and imposing the syncrhonous coordinate condition, 
a set of equations for the ordinary matter and the cosmic fluid perturbations is obtained: 

h + 2-h = -6-A m , (48) 
a a 

A m + (l + a)(tf-^) =0 , (49) 

2 

a Ti 

(l + a)^ + (l + a)(2-3a)-^-— aA m = , (50) 

a a 2 

A s + jj(e~) = o , (51) 

2 

CL Ti 

e + 3-e + — a s = o . (52) 

a 2cr 

The integration procedure is standard [Tj|, so we will just present the final results: 
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model 


perturbations 


material phase (£ — » 0) 


string phase (t — > 00) 


a = 






constant 


A s 


A (01X3(7177) + a 2 Is(n7])) 

rj3 2 2 




a = 1/3 




4-/ 7^ 


c\J\(nrj) + c 2 J_i(nrj) 


dr\ 


cosf-% hit) 


A s 


-^-(di/i (777/) + d 2 Ki(nr])) 

rfl 2 2 




a = 1 


A m 


1 f - 


CiJo(t7t/) + c 2 J (nr]) 


drj 


cos(nlnt) 


A, 


diI Q {nrj) + d 2 K (nr])) 





We must exhibit the behaviour of the perturbative model for t — > and t — > 00 to 
find an instability in the density perturbations of cosmic string fluid. Indeed, for small 
wavelengths, i.e. n — > 00, the modified Bessel functions I v (x) and K u (x) go as e ±x /^Jx. 
However, the ordinary matter has a very regular behaviour for t — > and t — >• 00. 

The instabilities presented above can be attributed to the matter content of the uni- 
verse and, mainly, to the approach used to describe the density perturbations. We will see 
that these instabilities do not appear in gravitational waves because the matter content 
of the universe plays a less important role, being only sensitive to the behaviour of scale 
factor. 

In fact, the same behaviour for the scale factor described above can be obtained by a 
scalar-tensor model, with a suitable potential, coupled to ordinary matter. For example, 
in the case of pressurelles ordinary matter (p — 0), the results for the background can be 
recovered if the potential reads 



V((f>) = V smh-\C(f)) 



Vo 



8ft 2 + 2 
a 2 



C 



V8TF+2 



n 



2vrGpo 



(53) 



where po is defined by the first integral of the conservation of the energy-moment tensor 
for the pressurelless fluid, p = We note that the scalar field representation for the 
two fluid model (pressurelles matter plus string fluid) requires a different potential with 
respect to the case of the one fluid model (only string fluid). 

At the perturbative level, the effect of replacing the barotropic fluid with negative 
pressure by a scalar-field model is quite the same as that described in section 3. With 
respect with the preceding table of solutions, the modified Bessel functions must be re- 
placed by ordinary Bessel functions for the perturbation in the exotic fluid. Hence the 
instability present in the small wavelength limit disappear again; the behaviour in the 
long wavelength limit is not changed. 



4.2 Gravitational waves 



Following the work [fL4}| , we have the perturbed equations that govern the evolution of 
gravitational waves 



h " - 2-ti + 
a 



fi 2 - 2 







a a 

where n 2 = (l/r 2 )(n 2 + 2k) and prime here denotes derivative with respect to 9 = rrj. 



(54) 
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After inserting the background solutions (^,^,^), equation ( |54|) can be generally 
rewritten in terms of a hypergeometric equation. Its final solutions for different phases of 
the evolution of the universe read as follows: 



model 



hi 



ho 



a 



l+x 



n -2+Vl-n 2 



X 



l+x 



-2-Vl-n 2 



X 



Ft 2- y/T 



h 2 - 

ll , 2 



n 2 



v ' l+x 



a = 1/3 



exp((vl — n 2 )r]) sinhr/ 



l + 2VT^, ^) 



exp(— (vi — n 2 )ri) sinhf? 



a = 



V 7 



iFi 



-1 - 
1 



l+x 



i+Va 1 



X 



ri 2 - 



v/4 



n 2 



7 
l 

2 



l+x 



X 



2 V / 4^ 2 , T ^ 



v 7 



2 



(-l+Vl-4n 2 )/2 



2 F 1 ^ + v / 4^,-l + v / 4 
1 + 2^4=7?^) 



n 2 



^1 



l+x 



(-l-v / l-4n 2 )/2 



X 



l-Vl-4n 2 l-Vl-in 2 



I-VT^W,^ 



w Z7- / l+Vl-4n 2 l+Vl-4n 2 
X 2^1 ( g ' 2 : 



where in the above expressions, 2F1 is a hypergeometric function and x = cos(r#), r being 
a constant. 

In this two-fluid model, the behaviour of gravitational waves in a closed universe 
exhibits, in what concerns the behaviour of the scale factor, the dynamic of an open 
universe. It would cause, also, distorsion in the spectrum of the CMBR. The determination 
of the evolution of perturbations and its connection with this distorsion is technically 
difficult to be evaluated. But, for the proposal of the present work, the fundamental 
feature to be retained in the above solutions is that they do not exhibit instabilities. 



5 Time dependent cosmological constant model 

We present for completeness the background equations based on the reference which 
is the traditional scalar-tensor theory with a potential that is equivalent to a time de- 
pendent cosmological constant model. This is one of the cases which can be represented 
phenomenologically by a fluid with an equation of state p = — | in what concerns the 
behaviour of the scale factor. The action is given by 

s = J V^[<PR - <tr x wsr%dv<l> + 2M(0)]^ + S ng . (55) 

We remark, however, that in the present case we have a non-minimally coupled scalar 
field, in opposition to the models described before. 
The field equations are 

g ^ rj~\ -j^ 

G^v = — -^ + A(0)^ + w0~ 2 (0, M 0,^--^0,a0, a ) + _1 (0;^-^n0) , (56) 
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^ + 2^ i - MM= 1 , T _^ K] _ (57) 



3 + 2u 3 + 2aA d<p 

We shall consider the case where oo = constant, A(0) = ap m and <p = <pit q , with c, m 
and g constants. The ansatz allow us to obtain analytical solutions in the form of power 
law, which leads to closed expressions for the perturbative equations. 

The energy-momentum tensor describes an ordinary perfect fluid, such that 

T^;„=0 . (58) 

The equations of motion are 



/ ,0* nA rn 8np LU(j) 2 d(j) 



-2*-^-!L + or = ^ + %f + U2*t. , (60) 
a a 1 a 1 <p 2 <p z <p a<p 

4> a<p _ 2c(l-m)0 m 87r(p-3p) 

<p a<p 3 + 2a; ~ 0(3 + 2w) 1 J 



The background solutions are 



model 


curvature 


a(t) 


<Kt) 


P 


A(f) 


p = 


ky^O 


ait 


-2 
(pit m 




A1A 2 


k = 


ait 


-'2 

0lt 1 ±v'3+2a, 




A1A 2 


p = ap 


any k 


ait 


lt -(l+3«) 


s / a 3(l+a) 


Ai/t 2 


p= -p 


any k 


ait 


0li 2 


p = const. 


Ai/t 2 



where ai, 0i, Ai, and s are constants. 

The solutions for the scale factor above are characteristic of an equation of state 
P<t> = — 3P0 • 

We remark that in all above solutions the scale-factor behaves as a oc t corresponding 
to a "coasting" universe. From the point of view of the background behaviour, the above 
solutions exhibit the same characteristics as the perfect fluid formulation. 

6 Perturbations in a time dependent cosmological 
constant model 

We extend our perturbative analysis for this model, first for density perturbations and 
then for gravitational waves. 
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6.1 Density perturbations 

The perturbed equations for the time dependent cosmological constant model are: 



a ■ 16n , . , . / 3au + 3a + u + 2" 

h + 2-h = — (A - A) )p 

a v ; V 3 + 2cu 



n 



+ ^A + 



2c(m 2 - l)0 r 
3 + 2cj 



■A = 



8tt 



(3 + 2u)(j) 



A(l-3a)p 



1 • 



A = (l+a)(-/ l -^, fc 



— (a 5 5f/ fe (l + a)p)= -a 3 apA.' 



(62) 



(63) 
(64) 
(65) 

where h = hkkO>~ 2 , A = S(p/(p and A = Sp/p. All functions are spatially expanded into 
spherical harmonics f(x,t) = f(t)Q, with V 2 Q = — n 2 Q. 

Next, we solve the above equations for the vacuum fluid (a = —1), the empty universe 
(p = 0), radiation phase (a = 1/3), and dust phase (a = 0). We use the background 
relations in order to simplify the final system of equations. 



6.1.1 Vacuum fluid phase (a = —1) 

In this case A = and we have the following equations: 

2 A A 

h+-h = (1 + 2cj)- + 8(1 + cj)- + 2A , 

6 6 6 



A + 7^ + 

6 



n n 2 lA 1; 
8 + — 



ai 



---h = , 

t 2 t 



with the solution 



A = t p and h oc t v , where p = — 2 ± y 8cj 
There is no divergent behaviour in the small wavelength limit. 
6.1.2 Empty universe (p = 0) 

Here, the system of second order differential equations is 



v 2 • 4 

h+-h = 

t m 



m + 2u + 2 



A_ _ 8(l + o;) A 
t 2 m t 



A + 



3m — 4 



m 



A + 



3 + 2u 

n 2 lln . 1 — m 
— -4 2 + m — 



+ 2A , 



A 1 /i n 

V m t 



(66) 
(67) 

(68) 



(69) 
(70) 
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whose solution is 



A = t p and h oc t v where p = — (l — m ± */9 + 6m — 3m 2 + 8c<j — m 2 ^—] , (71) 

m V V °i / 

where m is arbitrary for fc ^ and m = 1 ± V3 + 2cj for fc = 0. As in the preceding case, 
the solutions are stable. 



6.1.3 Radiation phase (a = 1/3) 

Combining conveniently the perturbed equations, we find a fifth order Euler's type equa- 
tion for A: 



\iv 

A- + ii + 

6 



4n 2 „„ „ , „ fc 
-— + 31-2rf-2- 
3 a x 2 af 



t 2 



16 n 2 ^ „ , „ fc 
— — + 22-8rf-8- 
3 ai 2 



+ 



1 4 

1 n n 
+ — 



3 ai 4 a x 2 V3 ' 3 



+ ld-2X)+2-4d-4 k 



a 2 



a 2 



A_ 



+ 



1 n 4 n 2 /2 fc 
3 a! 4 ai 2 V3 



a? 



A_ 



A_ 



= ,(72) 



with d = 3 + 2uo. Due to the residual coordinate transformation freedom characteristic of 
the synchronous coordinate condition, A oc t^ 1 is a solution of this equation. Hence, we 
can reduce it to a fourth order equation. Supposing a solution of the type A oc t r , we find 
the polynomial equation for r 



r 4 + 



whose solutions are 



4 n 2 „ , J\ , /In 4 



2d-2-)r 2 + 

3ai 2 af 



n 2 2 fc 
VW~^ ( 3 +2 ^ } 







(73) 



r ± 2 = 



2 n 2 

■-— + d+-± 

3 ai 2 af 



,2 fc 



9a! 4 v 3af 



^2 ^ 

2d)— + (d+ -o) 



a 2 ' 



(74) 



In the small wavelength limit, these solutions display an oscillatory behaviour, hence 
stability. In the longwavelength limit, this expression reduces to r± = ±w2(d+ ^). 



6.1.4 Dust Phase (a = 0) 

In this case, the equations governing the evolution of density perturbations are 
A + 2 j - (2 + uo)R^ = X-2(l + u)j-(2 + u;)S^ , 



A \n 2 



—; + S 

ai 2 



A A A 

T 2 ~ ~1 + R T 2 ' 



(75) 
(76) 



u,uu , - a \ +2uJ — • These equations may be solved supposing 
A = A t r and A = \ot r , r obeying the polynomial equation 



where R = a \ , » and S 



r 4 + 2r 3 + + S - (3 + u)R - (3 + 2cu))r 2 

/ n 2 \ n 2 

+ (—-(l+Lo)S+(3 + 2uj)R)r-(2 + uj)R— = , 
\Qi / ai 



(77) 
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This polynomial equation has no simple closed form solution. But, numerical integration 
reveals the stability of the model in the small wavelength limit. For example, fixing a\ = 1, 
choosing u — 1, k — and n = 1 we find the roots r\ ~ — 1.16±2.77i, r2 ± ~ 0.16±0.33z, 

while for n = 100, keeping unchanged the other parameters, we find r\ ~ 0.50±100i, 

r 2 ± 0.50 ± 0.97i. 



6.2 Gravitational waves 

In this case, the perturbed equation that govern the gravitational waves is 



,(f> a, 



n 



A 2 



^7+4- + 2 



w{a- 1) - 1 
3 + 2u 



a 2 a 2 2 <p 2 a<f) 



. 1 + m + ct>(l + a) 

- 2i sT^ "' 



h = Q . 



The solutions of the above equation are 



p = 
A 1 

a = — 1 

a = 
1 

a = - 



2 + m 



m 

h = C±t 
1 

3 + 2cj 



5 = -^(m 2 - m- 2- 2w) 



±,/-(^+A 2 ) 



3 - 8c; + 4cj 2 - 3— (1 + 2cj) 



3±. /9-4Aq-4- 



Am 



4A; 
ai 2 



2+ A k ., " 2 , 



(75 



(79) 

(80) 

(81) 
(82) 



where C± are integration constants. 

Here, the solution for the gravitational waves are also well-behaved and stable. 



6.3 The General Relativity limit 

In general, the Brans-Dicke theory reduces to the General Relativity theory when uj — > oo, 
except in some special cases, for example, when the trace of the momentum-energy tensor 
is zero [20, 21]]. The solutions described above do not have a well-behaved limit when 



lu — » oo. In fact, an inspection of the background equations shows that all solutions 
become trivial in that limit. We could expect that in this case, the Brans-Dicke field could 
become constant and the model reduces itself to General Relativity with a cosmological 
constant. But, the imposition that a oc t breaks this equivalence. 

Concerning the perturbed solutions, they become divergent when u — > oo. This only 
express the fact that the background scenarios have no sense in this limit. For finite u, 
the perturbed solutions exhibit growing and decreasing modes as usual, for scalar and 
tensorial perturbations. One important feature of the solutions found before is that when 
n — > 0, all dependence of the solutions on the wavelength of the perturbations is carried 
out by the integration constants, which must be fitted correctly in order to reproduce the 
power spectrum observed today. 
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7 Conclusions 



In spite of the fact that fluids of negative pressure have become crucial to understand 
many theoretical and observational aspects of modern cosmology, they are plagued with 
instabilities in the small wavelength limit. These instabilities appear mainly when the 
barotropic equation of state is such that the strong energy condition is violated. In this 
paper we have exploited the possibility that these instabilities are due to the hydrody- 
namical representation. This possibility was suggested by the fact that, while density 
perturbations are bad behaved in the small wavelength limit, gravitational waves are well 
behaved in the same limit. Since gravitational waves depend closely on the behaviour 
of the scale factor but are quite insensitive to the matter representation, the instabilities 
should be due to the fact that in the hydrodynamical representation, the equation of state 
is fixed for all wavelength, while in a more fundamental approach it could depend on the 
scale of the perturbation. 

As a matter of fact, this possibility was first suggested in |1C|, and it has been studied 
under certains conditions in (T^|. In reference |T(J, this problem has been briefly treated 



in the realm of minimal scalar-tensor model which was intended to cope with the dark 
matter problem through the introduction of a scalar field with a convenient potential. 
However, if we are interested in a field that can drive the accelerated expansion of the 
Universe, as it is the case in this work, the energy of the scalar field should have a smooth 
distribution, since it should not be present in the local clusters. Such smooth distribution 
can be obtained considering that the pressure associated to this field is negative, such 
that it does not clumpsy in large scale; but in order to avoid instabilities at small scales, 
the effective equation of state should become positive in this limit, and small fluctuations 
in this field should oscillate like an accoustic wave. 

In the present work, we have verified that, when a fluid of positive pressure can mimic 
a scalar-tensor model, both formulations are essentially similar in the background and 
perturbative level. However, when we consider a fluid of negative pressure, the equivalence 
exists only at the background level: at perturbative level, the models behave in a complete 
different way. In particular, there is agreement between both representations only in the 
large wavelength limit: for small wavelengths, the hydrodynamical model is unstable, 
while the corresponding scalar field representation exhibits accoustic oscillations. Hence, 
in situations where negative pressures are concerned, a field representation leads to a 
much more complete scenario, being closer to a realistic model. 

This can be understood by remembering that when we pass from a hydrodynami- 
cal representation to a fundamental one based for example on scalar fields, we change a 
relativistic Euler's type equation to a Klein-Gordon equation: the sign of the laplacian 
operator in these equations are the same only when the hydrodynamical pressure is posi- 
tive; otherwise, it changes sign and, in the perfect fluid model, accoustic modes give place 
to exponential modes, resulting in the presence of instabilities. In the long wavelength 
limit the laplacian operator can be neglected and the results agree in both representation. 

We must notice, however, that the scalar-tensor model that corresponds to a given 
perfect fluid results is quite model dependent. For example, the one fluid model gives 
a potential that is different from that of a two-fluid models. It would be interesting 
to employ in a two fluid model, which is closer to a realistic situation, the potential 
obtained in the one fluid model, in spite of the great complexity of the equations. We can 
speculate if the resulting effective equation of state evolves in a quite similar way as the 



16 



usual quintessence models. 

We have extended this study for the case where the cosmic string fluid mimic a variable 
cosmological constant model, derived from a non-minimal coupled scalar field with a 
potential. The conclusions are basically the same as before. But, we must stress that, 
in this case, there is a Brans-Dicke type coupling parameter u that does not recover the 
known General Relativity solutions when uj —>■ oo. Moreover, the perturbative behaviour 
may become unstable in this limit, for any scale. But, in general, for finite values of u the 
perturbations do not exhibit either the instabilities that are present in the corresponding 
hydrodynamical model. 

Finally, the fact that both approaches give the same behaviour for the long wavelength 
limit even if the pressure is negative, implies that the predictions for the power spectrum 
of the anisotropy of the CMBR for small values of the multipolar expansion parameter I, 
that is, for very large structures, is not spoiled by the employment of the hydrodynamical 
representation. However, for large values of I, that is, small angular separation, the 
employment of a field representation, mainly when negative pressures are involved, seems 
crucial. 
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